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Abstract 

By means of Ito calculus it is possible to find, in a straight-forward way, the 
analytical solution to some equations related to the passive tracer transport problem 
in a velocity field that obeys the multidimensional Burgers equation and to a simple 
model of reactive tracer motion. 
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1 Introduction 



In a recent paper JlJ Saichev and Woyczynski have obtained exact solutions in arbitrary 
dimensions of equations of hydrodynamic type related to the Burgers equation (see, e.g., 
@, 0, [|J, ||) for an irrotational velocity field with two models of coupled passive or 
reactive tracers. For the inhomogeneous Burgers equation and more general models of 
tracers, closed but less explicit solutions are obtained as path integrals. 

The main idea of is to reduce the system specified by the forced Burgers equation 
together with the advection-diffusion-reaction equation, to a pair of coupled linear diffu- 
sion equations with variable coefficient which can be analytically solved by means of the 
Feynman-Kac formula (see, e.g. 0). 

They also show that the same methodology can be used to generate exact solutions 
of a non linear reaction diffusion equation coupled with a Burgers-like velocity field also 
depending on the concentration. 

The Feynman-Kac equation expresses the solution of a parabolic PDE without drift 
in terms of a conditional average over Brownian trajectories. A generalization of the 
Feynman-Kac formula when a drift is present, is supplied by the Cameron-Martin-Gir- 
sanov formula ||. This observation allows to recover the results in |J without using any 
auxiliary field. 



2 The Girsanov formula 

In this section a generalized version of the Girsanov formula is recalled and its relevance 
to parabolic PDEs is explained. 

Theorem 1 Let x^\ x^ E H d be solutions on the interval < t <T of the stochastic 
differential equations 

dx^it) = P\x®(t),t)dt + a(x®(t),t)dw(t) 
fW(o) = x i = l,2 (1) 

where w(t) is a d-dimensional Brownian motion, b®(x,t) (i = 1,2) and a(t,x) are re- 
spectively Borel-measurable, R d -valued functions on [0,oo[xR d and a Borel-measurable, 
[0,oo[xR d x R d -valued function with bounded inverse Vx,t. If b^'(x,t) (i = 1,2) and 
a(t,x) satisfy the assumptions of the existence and uniqueness theorem for the solutions 
of (W, then the probability measure fi 2 of x^ will be absolutely continuous w.r.t Hi of x^ 
and 

$%«(*)) = e A > 
d/ii 

At = f a(x w (s),s)-dw(s)-- t \\d(x {l \s),s)\\ 2 ds (2) 
Jo 2 Jo 



where 



2 Jo 

a(x, t) = a' 1 (x, t) [b {2) {x, t) - (x, t) } (3) 
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see H pag. 279 for proof and details. 

A straight-forward consequence is that for any reasonably smooth function fk{x), with 
k ranging from 1 to n and for all n-tuples (tx, ■ ■ ■ , t n ) such that < tx ■ ■ ■ < t n < t we 
have: 

n n 

E*{Ufk@ W (tk))} = E 3 {U h{x {1 \t h ))e A *} (4) 

fe=l k=l 

This result can be exploited in order to write the solutions of parabolic PDEs as path 
integrals on Wiener trajectories (see e.g. 0). 

In the following sections the result @) will be applied, disregarding the conditions 
on the drift field, in order to derive formally the announced results. The advantage of 
this approach is to supply a direct physical interpretation for the solutions in terms of 
stochastic trajectories. 



3 Simple applications 



As a first application let as consider the d-dimensional homogeneous Burgers equation 
with rotation free initial condition: 

d t v + v ■ Vv = uAv 

v(x,0) = V9 (f) (5) 

The physical meaning of ([5|) is that the velocity field is, on the average, constant along 
the trajectories generated by the stochastic differential equation: 



dx(s) = —\7<d(x(s),t — s)ds + \ / 2vdw(s) 
f(0) = x (6) 

If we introduce: 

dz(s) = \plv dw(s) 
z(0) = x 

z(t) l = Af(x,2vt) (7) 
we can exploit Girsanov's theorem and write: 

V9(f,i) = E s {V z{t) Q,{z{t))z~ Zt } (8) 

where: 

Z t = ^= [ t VQ(z(s),t-s)-dw(s) + — f \\Ve(z(s),t- s)\\ 2 ds (9) 
V2 v Jo Av Jo 

We can eliminate the stochastic integral in (|9|) by means of 
d s Q(z(s), t — s) = 

= {-d t _ s G(z(s),t- s) + p Ae(z(s),t - s)} ds + V2uVe{z(s),t - s) ■ dw(s) (10) 
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This substitution is useful since the potential satisfies: 

e(x,o) = e (f) (n) 

So we get to 

Ve(f,t)exp[-^M] = ^{V z(t) e (^))exp[-^S^]} {12 ) 

Finally, we can derive the explicit expression for the velocity potential by means of a 
simple integration by parts and by exploiting the homogeneity of ([?[) : 

Q{x , t) = -2v In E s {eM~ }} (13) 



which means: 



where 



oo 



e(f,t)=-2,M/ e-^ ¥ ^ w ] (14) 



Hx,y,t)= { -^^ + e (y) (15) 



If we consider a class of initial conditions such that 



\x\\ 2 



for ||x|| | oo (16) 



then (|TJ|) is always well defined and, in the limit t j 0, it is consistent with the initial 
condition 0(x, 0) = Oo(^)- Therefore, we have recovered the already known result of the 
Hopf-Cole theory (see e.g. H, M) stating that the solution of the Burgers equation with 
irrotational initial condition is given in any dimension at arbitrary time t by 

v(x,t) = V0(f,t) (17) 

Let us now consider the system: 

d t C + v-VC = fiAC + VC + g 

C(f,0) = C (x) (18) 

where v is given by (|5|). The external drift V and the volume force g are functions of x 
and t; V, g and the initial data Co are smooth functions growing, as ||x|| goes to infinity, 
more slowly than ||x|| 2 . In the general case this equation can be formally integrated as a 
path integral in the form (see e.g. 0): 

C(x,t) =E s {C (z(t))eM-Zt+ f t V(z(s),t-s)ds]}+ 

J o 

+ E s {[ g(z{s),t- s)exp[-Z s + [' V(z(u),t - u) du]ds} (19) 
Jo Jo 



where z(t) and Z t are given respectively by (0) and @. 

If fi = v and V = g = the situation is more simple and the solution can be easily 
expressed as an ordinary integral. 



C(x,t) =E x {C (z(t))e- Zt } (20) 

It is now sufficient to proceed, as before, to the elimination of the stochastic integral in 
(§) by means of fllPP and to use flllD to arrive to 

r 6(x,t) f 00 ${x,y,t) d d y 
C(x,t)=eM-^ r ]J_ oo C (y)eM — ] ^— ^ (21) 

or, in more explicit terms: 

which is the first result of Saichev and Woyczynski. 

If the diffusion coefficient in (|18|) is fi ^ v the velocity and the concentration fields 
are, on the average, constant along different stochastic trajectories, so we can no more 
use ([0]) to reduce the path integral to a finite dimensional one. Still, we can formally 
write: 

C(x,t) = e^^{C (r(t))exp[- e ° (z/(t)) - f A9(i», t - s) ds]} 

I [i 2 (i Jo 

z\t) l = Af(x,2fit) (23) 



where 9 (x, t) is given by (|Tj) . The expression (|23|) holds true as far it is well defined (non 
divergent). 

4 A reaction-diffusion model 

Another possible simple application of the Girsanov theorem is the solution of the follow- 
ing reaction diffusion model considered in 

d t v + v-Vv = vAv + 2vkVC v{x , 0) = V9 (f) (24) 

d t C + v ■ VC = is AC + kC 2 C{x,0) = C {x) (25) 

where k is a constant. 

Let us consider the stochastic differential 

ds { e kJoC(x(u),t-u)du C ^ s y f _ s )} = 

= e k fo c ^ t -^ du {kC 2 (x(s),t-s)ds + d s C(x(s),t-s)} (26) 
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where 

dx(s) = -Ve(x(s),t -s)ds + V2vdw(s) (27) 

If we take the expectation value of (|6|) and equation (fj4|) is satisfied, then it is easy 
to see that: 

C(x, t) = E s {C (z(t)) exp[-Z t + k [ C{z{s),t - s) ds}} (28) 

Jo 

where z(t) and Z t are given respectively by (^) and (§) 

Again we can use (|T0| ) and we obtain an expression of the form fl2~T|) where, now, G(x, t) 
is given by solving ([24]). This can be done by observing that, in terms of the potential, it 
has, according to the notation in the form: 

d t Q + -Ve- VQ = vAG + 2vke& a (29) 



Here we have 



f(s,y,t) d d y 



«.<*.«) =y_ £*.«>«- - < 3 °) 

and $(x, y, t) is defined by (|i~5|). 

The Ito calculus suggests that equation ( p9|) can be rewritten by means of the stochas- 
tic process (|7|) as 

ff g Rexp[- e(z(s) '*~ 5) ]} = kE*{a (z(s),t-s)} (31) 

which implies 

Q{x,t) = -2u lnE^ie- 21 ^ - kta o (z(t),0) + C s d s a (z(s),t - s) ds} (32) 

Jo 

Since ao(x, t) is solution of the heat equation we have 

E?{[ sd s a (z(s),t- s)ds} = (33) 
Jo 

Therefore, we can conclude that the velocity potential is now 

Q(x , t) = -2u ln[ b (x, t)-kt a (x, t)] (34) 

with 



roo <s,(3,y, t ) d d y 
b (x,t)= I e *» {A7riyt)d /2 ( 35 ) 



— oo 



while the reactive tracer is 



Oq(x, t) — kt ao{x, t) 



which is the second result of Saichev and Woyczynski. 
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5 Conclusion 



By means of Ito calculus is it possible to obtain the result reported by Saichev and 
Woyczynski in (PJ) in a natural and straight-forward way. 
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